Weak Wave Turbulence Scaling Theory for Diffusion and Relative Diffusion in 

Turbulent Surface Waves 
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We examine the applicability of the weak wave turbulence theory in explaining experimental 
scaling results obtained for the diffusion and relative diffusion of particles moving on turbulent 
surface waves. For capillary waves our theoretical results are shown to be in good agreement with 
experimental results, where a distinct crossover in diffusive behavior is observed at the driving 
frequency. For gravity waves our results are discussed in the light of ocean wave studies. 



I. INTRODUCTION 

The study of particles moving on surface waves has 
shown that the particle motion often is far from being 
Brownian [jl], ||, [| ||, |j| . Similar results are found in ocean 
studies The turbulence observed in surface waves is 
strongly influenced by the presence of a dispersion rela- 
tion (for deep-water surface waves) lo\= gk + ak 3 /p be- 
tween the wave oscillation frequency and the wave vector 
k . Here g is the gravitational acceleration, a is the sur- 
face tension coefficient, and p is the density. Experimen- 
tally surface waves can be studied by vertically oscillating 
a fluid with a free surface. On the surface of the fluid, 
waves with a clearly discernible wavelength A are then 
formed if the amplitude of oscillations exceeds a critical 
value, the so-called Faraday instability ||. The excita- 
tion of surface waves in the Faraday system is a paramet- 
ric effect giving rise to a fundamental wave frequency 
equal half the driving frequency. The frequency Q and 
wavelength A are associated through the dispersion rela- 
tion. At small frequencies, k < k a — {gp/a) 1 / 2 , the first 
term in the dispersion relation dominates, and the waves 
formed are called gravity waves (typical of ocean waves). 
At large frequencies (k > k a ) the effect of gravity can be 
neglected, and the second term dominates. The waves 
thus formed are called capillary waves. For water, the 
wave length a = 2ir/k a , separating gravity and capillary 
waves, is about 2 cm. 

Several studies have been carried out to measure the 
diffusivity of particles moving on capillary surface waves 
generated above the Faraday instability |l|, ||, [| [| [|, 
|9j. Two particles separated by a distance R move more 
rapidly apart as R is increased. The relative diffusivity 
can be defined as 



, dR 2 , 



(1) 



where Sv is the longitudinal component of the velocity 
difference. The variance er 2 of the relative diffusivity 
dR 2 /dt is related to the velocity difference squared, 



r) 



(\Sv(R)\ 2 ) 



(2) 



The purpose of the present investigation is to exam- 
ine the applicability and limitations of the weak wave 
turbulence theory in explaining experimental scaling re- 
sults obtained for the diffusion and relative diffusion of 
particles moving on turbulent surface waves. Diffusion 
theories based on random wave fields have likewise been 
considered @, and Herterich and Hasselmann Jll[ par- 
ticularly considered diffusion by surface gravity waves. 
For relative diffusion on capillary waves we find the ex- 
ponent — 1/4 as in experiments in the regime R > A. 
For gravity waves, we find an exponent 7/6 ~ 1.14, a 
result which is intriguingly close to the value 1.15 found 
in ocean studies by Okubo Q . We note that the relative 
diffusion exponent 4/3 obtained in fully developed turbu- 
lence theory previously have been used in the discussion 
of experimental data |l^] . 

To obtain the diffusivity, the particle displacement 
Ax(t) = x(to + 1) — x(to) over a time t is measured along 
an arbitrary axis for many initiation times t$. From the 
resulting distribution the variance V(t) — ([Ax(t)] 2 ) is 
found, and the diffusivity D can be extracted, V{t) = 
2Dt. For Brownian motion the diffusivity D is a con- 
stant. For deep water capillary waves (we only consider 
deep-water surface waves) formed at high frequencies, the 
diffusion has been found experimentally to be anomalous 



V(t) ~ t x , 



(3) 



It has been found experimentally that in capillary waves 
(f ) 2 oc 1 - aOR/A)- 1 / 4 , with a = 0.66 @. 



with A > 1 (super-diffusive) || [|, ||, g|. The exponent is 
observed to change drastically from a value A = 1.6 — 1.9 
at length scales below A to a value A = 1.0 — 1.3 at 
length scales above A. The values of A is observed to 
decrease with increasing drive; at large drive A attains 
the Brownian motion value [|. Assuming a connec- 
tion between spatial and temporal correlations we find 
for capillary waves A = 16/9 for length scales below A, 
and A = 1 for length scales above, i.e. values rather close 
to the experimental values. Ocean studies have been per- 
formed using floaters [^3), chemical tracers Q, and near 
surface drifters. To what extent these studies relate to 
surface waves is unclear. Our results do, however, indi- 
cate a connection at least for the chemical tracers, where 
Okubo finds an exponent of 2.34 and from weak wave 
turbulence theory we find 8/3. 
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II. SPATIAL CORRELATIONS AND RELATIVE 
DIFFUSION 

Spatial correlations between pairs of particles at r and 
r + R are related to the velocity difference squared 

(\5v(R)\ 2 ) = <|v(r + R)-v(r)| 2 ) 

= 2[(|v(r)| 2 )-(v(r + R).v(r)}] . (4) 

Following the calculations of Ref. fL4|| , we have an ex- 
pression for the space-correlation of the velocity 



C(R) = (v(r)v(r+R)) 



kui k n k e 



ikR cos ( 



dk 



k 2 u) k n k J (kR)dk ,(5) 



where n k is the isotropic Kolmogorov spectrum 15 , and 
Jo is the zeroth order Bessel function. 

For capillary waves, the dispersion relation is tu k — 
(a I 'p) x / 2 k a with a = 3/2, and the isotropic Kolmogorov 
spectrum (three wave interactions) is |l5, ttq| 



pl/2 3/ V -l/4 fc - 



(6) 



where — 17/4 and P is the energy flux assumed to be 
constant. From measurements of the spectrum of wave- 
amplitudes (jl7| [l8| we know that the spectrum levels off 
below K = 2tt/A. In order to account for this change 
we introduce a scaling function g(k/K) across the drive 
multiplying n k . At k < K, g(k/K) must rise from zero 
in such a way as to quench the decrease of n k . At k > K , 
g(k/K) can be assumed constant (up to the inverse dis- 
sipative scale kd where the theory breaks down). The 
motion at length scales below the dissipative scale is ex- 
pected to become ballistic. Also the system size L is a rel- 
evant length scale. Obviously the diffusion change char- 
acter near the system size since particles cannot travel 
longer distances than L. Also we shall consider the char- 
acteristic length scale a separating capillary from gravity 
waves. 

The integral for the spatial correlations is given by 



C(R) 



•„2-t-a-/3 



g(k/k F )J Q (kR)dk 



(7) 



where the specific k dependence for io k and n k has been 
inserted. For capillary waves 2 + a — (3 = —3/4. This 
integral depends explicitly on R. Consider the value 
Kr ~ 2.4/i?, the lowest k value at which Jo(kR) = 0. 
For Kr > K, the dominant contribution to the integral 
comes from the range K < k < Kr, where the scaling 
function g is constant, 



C{R) 



k 2+a -Pg(k/K)J a {kR)dk 



u3+a-0\K R 

K \K 



(8) 



The result for the velocity difference is 
(\Sv(R)\ 2 ) ~ 1 - b(KR)- 3 - a+fi ~ 1 - b(KR)- 1 ^ 4 , (9) 



where & is a constant of order one [Tg|]). 

For Kr < K, the Bessel function oscillations set in 
below K while the integrand is still increasing in size due 
to the behavior of the scaling function g{k/K). Therefore 
the dominant contribution to the integral comes from a 
peak centered around K, it is oscillatory in nature with 
a vanishing envelope falling off like R?/ 2 . For large R 
values we essentially have (\Sv(R)\ 2 ) constant. 

For gravity waves the dispersion relation is uj = g 1 / 2 ^ 
with a = 1/2 (four- wave interaction), and we have energy 
flow towards higher frequencies compared to Q giving rise 
to the isotropic Kolmogorov spectrum |15L B(| 



(10) 



where (3 = 4 and P is the energy flux assumed to be con- 
stant. Below !7 a constant wave- number flux Q towards 
lower frequencies yields [pTI 



(11) 



with /3- a/3 = 23/6. 

As above we invoke a scaling function g(k/K) describ- 
ing the change in the behavior of n k as k crosses the 
wavenumber K of the drive. According to the form of 
the spectra we have 



g(k/kj 



1 if K < k < k a , k d 
k a / 3 iik L <k<K. 



(12) 



Other corrections are relevant at length scales outside 
the above regime. The motion at sufficiently large length 
scales (fc < ki,) may be expected to become Brownian, 
as also found in studies of motion of drifters near the 
surface of the oceans 0. 

For the velocity difference squared we have 



(\6v(R)\ 2 ) 



k 2 +<*-Pg(k/K)[l - J Q (kR)]dk 

k~ 3 / 2 g(k/K)dk . (13) 



k r 



For Kji > K, the scaling function is constant in the 
integration region above K^, and we have 



(IM^)I 2 )-^ 



-1/2 



(KR) 



1/2 



(14) 



For Kr < K the part of the integration, where g(k/K) 
fc 1 / 6 , gives the R dependent contribution, and 



(\Sv(R)\ 2 ) + const ~ Kj 



-1/3 



(KR) 



1/3 



(15) 



The growth of (\Sv(R) \ 2 ) stops at length scale L. We note 
that the exponent 1/3 for the velocity difference squared 
implies an exponent 7/6 for the relative diffusivity mea- 
sure a. It is interesting that Okubo H finds the rela- 
tive diffusion exponent to be 1.15 for oceanic diffusion, 
although we should emphasize that the relative diffusion 
exponent 4/3 obtained in fully developed turbulence the- 
ory has generally been applied to explain experimental 
data 111. 
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III. TEMPORAL CORRELATIONS AND 
CAPILLARY WAVES 

Next we consider single particle diffusion in capillary 
waves. For arbitrary times, the mean-square displace- 
ment (r 2 ) = (\r(to + t) — r(to)\ 2 ) is expressed in terms 
of the Lagrangian velocity correlation function C(r) = 
(v(r(t))v(r(t + r))> as gg 



(r 2 ) =2t (1 - T/t)C{ T )dT 



(16) 



In the Taylor limit of very small t, C(t) ~ (v 2 ), where 
(v\) is the mean-square Lagrangian velocity, to be deter- 
mined by sampling along particle orbits. In this case 



(r 2 ) « (vl)f 



(17) 



For incompressible flows, the mean-square velocity ob- 
tained by Eulerian sampling, {v 2 E ), is identical to the 
mean-square velocity, (w 2 ), obtained from Lagrangian 
measurements |23|. The subscripts on (v 2 ) will be omit- 
ted from here on. For large times, t ^> tl, where 
t l = j^zj Jq 00 C(t)At is the Lagrangian integral time 
scale, Eq. ([l6[) gives a diffusion-like dispersion 



(r 2 ) ~ 2(v 2 )r L t 



(18) 



First we consider capillary waves. If we invoke a scaling 
form (r 2 ) ~ t x at t < T, then Eq. (Ill) gives 



C(r) 



A-2 



(19) 



On the other hand, invoking the same scaling (r 2 ) ~ t x 
and using Eq. d|) yields 



C(t)~(t*)-* 



(20) 



The result is a self-consistency condition for the exponent 
A. We find 



A 



5 + a - (3 



16 



(21) 



This value obtained for A is in good agreement with ex- 
perimental observations which vary ±10% from this value 
H [Ij. From this the exponent for the time correlation 
function becomes —2/9. 

Above it has tacitly been assumed that the time scale t 
is larger than the one corresponding to the dissipative fre- 
quency cod corresponding to the inverse dissipative scale 
kd- Experimentally, the sampled diffusion distances at 
small times may approach the dissipative scale at low 
values of the drive. Below this scale A has the ballistic 
value 2, and the observed slope (1.8-1.9) slightly larger 
than the theoretical value 16/9 may be an effect hereof. 

We have also assumed that the particles can follow 
the fluid velocity field, i.e. if v is the velocity field of 
the particles and u is that of the fluid, our assump- 
tion is that v = u, or at least that the weak ansatz 



C(r) = (u(t + r)u(i)) is fulfilled. However, at large driv- 
ing amplitudes the particles may not follow the fluid ve- 
locity. As an example, assume that corrections for the 
velocity correlation takes the following form 



(v(i)v(t + r)> 



(U(i)u(t + T)) 



(22) 



where v and u are the velocity fields of the particles and 
of the flow respectively. Following the calculation above, 
we then find 



A 



4-25 
5 + a - (3 



16-8<S 
9 



(23) 



which may explain why a slightly smaller exponent than 
16/9 is observed for higher driving amplitudes (e.g. S = 
0.2 gives an exponent A = 1.6). 

Equivalcntly, the self-consistent condition for A 
(Eq. (pl|)), can be obtained in terms of the energy spec- 
trum 



E k ~ k~ v , 
where the energy per area is 



E k dk 



E k dk 



kuknkdk 



we have the relationship v = — (1 -\-a — 0), and 
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A = 



(24) 



(25) 



(26) 



The capillary wave exponent is v = 7/4. 

For t > T, the time correlation C(r) is oscillatory and 
therefore negligible for t > T. In this case, Eq. ( |lS| ) 
applies 



(r 2 ) - t , 



(27) 



and A = 1. The Brownian motion value obtained for 
A at times r above T is in agreement with experimental 
observations for large drives. The slightly larger values of 
A at lower drives may arise from corrections to scaling. 
Experimentally, at sufficiently large times (above those 
considered in |], [|), the sampled diffusion distances are 
limited by the systems size L and eventually one must 
have A = 0. 



IV. TEMPORAL CORRELATIONS AND 
GRAVITY WAVES 

For gravity waves, we can follow a similar assumption 
as the one proposed for capillary waves. However, in 
contrast to the situation for capillary waves, the scaling 
is given for (lv(t) — v(t + r)| 2 ). Thus, for r < T, we have 
from Eq. (Q 



(|v(t)-v(t + r)| 2 ) 



-A/4 



(28) 
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and for r > T, we get from Eq. (|15|) 

<|v(t)-v(t + r)| 2 )^rV6. 

Assuming 

i? 2 

the exponents can now be calculated using Eq. 
the self-consistent condition 



,A-1 



i 5A / 8 small t 
t 7A/i2 large t 



(29) 

(30) 
from 

(31) 



The values obtained for A are 8/3 for small t and 12/5 
for large t. Interestingly, a similar situation is observed 
in ocean studies where A = 2.34 is found in several dye 
studies collected and analyzed by Okubo K9. 



V. CONCLUSIONS 

Good agreement is found between theoretical and ex- 
perimental scaling results for the diffusion and relative 
diffusion of particles on weakly turbulent surface waves. 
In the case of relative diffusion, for capillary waves the 
exponent -1/4 for the space correlation function is recov- 
ered. For gravity waves, a scaling behavior is derived for 
the velocity difference squared with the exponent chang- 
ing from 1/2 at small distances to 1/3 at large distances. 



This corresponds respectively to scaling exponents 5/4 
and 7/6 for the relative diffusivity. Intriguingly, Okubo 
finds a relative diffusion exponent of 1.15 for oceanic dif- 
fusion ||. This differs from the relative diffusion expo- 
nent 4/3 obtained in fully developed turbulence theory, 
previously used in the discussion of experimental data 

& 

For capillary waves the theory yields a crossover in the 
diffusivity at the wavelength A with a change in diffusion 
exponent A from A = 16/9 ~ 1.78 to A = 1, in good agree- 
ment with experimental findings. Particles suspended 
on the fluid surface may not follow the fluid far above 
the Faraday instability, and we suggest that this may 
explain the small discrepancy between weak turbulence 
theory and experimental results in this regime. For grav- 
ity waves the weak turbulence theory gives a diffusion 
exponent of 12/5 = 2.4 for long times which compares 
favoritely with the exponent 2.34 found by Okubo for 
dye diffusion. To what extent the weak wave turbulence 
theory really explains the oceanic diffusivity is unclear 
but it is fascinating to speculate that wave motion tak- 
ing place mostly on a short scale could be responsible 
for the long time diffusion in oceans. At very long time 
scales, i.e. above the Lagrangian time scale, ordinary 
Brownian motion is expected. 
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